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Abstract 
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1 Introduction 

\ We address the problem of the separation of variables for the Hamilton- Jacobi 
equation within the theoretical scheme of bi-hamiltonian geometry. The main 
aim is the construction of variables of separation for the given integrable system 
without any additional information (Lax matrices, r-matrices, links with soliton 
equations etc.) 

The paper is organized as follows. In Section 2 we determine the Bogoy- 
avlensky system on so(4). In Section 3, the necessary aspects of bi-hamiltonian 
geometry are briefly reviewed. Then, we discuss a possible application of these 
methods to calculation of the polynomial bi-hamiltonian structures for the given 
Bogoyavlensky system. In Section 4, the problem of finding variables of separa- 
tion and corresponding separation relations is treated and solved. 

All the computations for the paper have been done by the computer algebra 
system Maple. It allows us to solve overdetermined polynomial differential and 
algebraic systems of equations. A major concern for any operations performed by 
the program is the complexity of resulting expressions. So, we can say that this 
note belongs mostly to so called "experimental" or "computational" mathematical 
physics. 
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2 The Bogoyavlensky system on so (4) 

Let 7 = (71,72,73) and M = (M 1; M 2 ,M 3 ) be the two vectors of coordinates 
and momenta, respectively. We postulate the following Poisson bracket on this 
six- dimensional phase space 

{M i5 Mj} = e ijk M k , {M h 7^} = e ijk 7*, {7^ 7^-} = x 2 e ijk M k . (2.1) 

Here Eij k is the totally skew-symmetric and x is a parameter. It is well-known 
that any linear Poisson bracket is defined by an appropriate Lie algebra. The 
cases x = 0, x 2 > and x 2 < correspond to the Lie algebras e(3), so(4) and 
so(3,l). 

The Poisson bracket (12.11) has the two Casimir functions 

d = M l7l + M 272 + M373 , C 2 = x 2 {Ml + Ml + M 3 2 ) + 7l 2 + 72 2 + 73 2 . (2.2) 

Hence, for the Liouville integrability of the equations of motion only one ad- 
ditional integral functionally independent of the Hamiltonian and the Casimir 
functions is necessary. 

The nontrivial class of quadratic homogeneous Hamiltonians of the form 

H = (M, AM) + (M, B 7 ) + (7, C 7 ), (2.3) 

where A, B and C are constant 3 x 3-matrices, has many important applications 
in the rigid body dynamics [3]. 

There are two classical integrable cases, one found by Chaplygin and one by 
Goryachev, where the additional integral of motion is of fourth degree. Namely, 
at x = and C\ = integrals of motion are in the involution with respect to the 
Poisson bracket (12. ip 

Hi = Ml + Ml + 2M 3 2 + ci( 7l 2 - 7 2 ) Cl G R, (2.4) 
H 2 = (M 2 -M 2 + ci7|) 2 + 4M 2 M 2 . 

It is so-called Chaplygin system on the sphere pE]. The Lax matrices and r- 
matrix formalism for this system have been obtained in [61 [8] , the corresponding 
bi-hamiltonian geometry has been studied in [121 US US] • 

Bogoyavlensky found pull-back of the Chaplygin system to the so(4) algebra 

H x = (l-x 2 ai)M 2 + (l-x 2 a 2 )M 2 2 + (2-x 2 a 2 -x 2 ai)M 2 + (ai-a 2 )(7 2 - 7 2 ), 

(2.5) 

H 2 = ((1 - x 2 a 2 )M 2 - (1 - x 2 ai )Ml - (a 2 - a^jfj + 4(1 - x 2 a 2 )(l - x 2 a 1 )M 2 M 2 . 

and proved that equations of motion can be integrated by means of elliptic func- 
tions [I]. 
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Remark 1 In fact we can found two integrable at G\ — systems on so(4) in 
the Bogoyavlensky book pp. The Hamilton function for the first system 



#1 = biMl + b 2 M\ + 6 3 M 3 2 
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x* 



J 3 2 

-7i 



2&i - 63 2 6 2 - ^3 + 6l 



X^ 



■72 + 



x^ 



73 (2-6) 



coincides with the Hamiltonian i?i (12. 5 p up to Casimir function 



if we put 

61 - 63 = x 2 a 2 - 1, 6 2 - &3 
The second system with the Hamilton function 

1 



x di — 1 



Hi = c x M{ + c 2 M 2 z + c 3 M 3 z + 



2x 2 



(c 2 + c 3 ) 7 r + (ci + c 3 )7 2 2 + (ci + c 2 )7 3 2 (2-7) 



has the following additional integrals of motion 



if,- 



\ 2 

(cj - c fc )7? + (cj - Ci)jl + (cj - c fc )7 4 2 J + 4(cj - Cj)(cj - c fc )7*7*, 



where {r/'fc} is one of the possible cyclic permutations of subscripts {123}. It is 
easy to prove that 

aKi + f3K 2 + 7^3 = 0, iff a + /3 + 7 = 0. 

According to [H] , Hamiltonians (12.61 ) and (12. 7\i are related by the Poisson map 

Mi -> — , 71 ->> xM l5 M 2 ->• i 2 -, 72 ->• xM 2 , M 3 -> M 3 , 73 -> 73, 

X X 



and change of parameters 

d = (262-63), c 2 = (26 1 -6 3 ), 



c 3 



^3 • 



3 The bi-hamiltonian structure 

In order to get variables of separation according to general usage of bi-hamiltonian 
geometry firstly we have to calculate the bi-hamiltonian structure for the given 
integrable system with integrals of motion H\ } 2 (12. 5 p on the Poisson manifold so(4) 
with the kinematic Poisson bivector P and the Casimir functions C x ^ ( 12.21) : 
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Following to [TU Q21 El [EE1 US] we are looking for solution P' of the equations 

{H x , H 2 }' = (P'dHi, dH 2 ) = 0, [P, P'} = [P', P'] = 0, (3.2) 

where [., .] means the Schouten bracket. 

Obviously enough, in their full generality equations (13. 2 p are too difficult to 
be solved because it has infinitely many solutions [TCJj [TI] . In order to get some 
particular solutions we will use the additional assumption 

P'dC h2 = 0, (3.3) 

and polynomial in momenta M ansdtze for the components P^ of the desired 
Poisson bivector P' = Yl Pij 9% A dj. 

Substituting polynomial ansdtze into the equations (13.2113.31) and demanding 
that all the coefficients at powers of M vanish one gets the over determined system 
of algebro-differential equations on functions of x which can be easily solved in 
the modern computer algebra systems. The computation was performed using 
the computer algebra system Maple and, in contrast with the Chaplygin case 
[To] , these calculations were performed non-automatically with essential manual 
interaction. In this way we get a lot of real and complex solutions, which will be 
classified and studied at a future date. 

In this note we will not consider a complete classification and restrict ourselves 
by discussion of one example only. In order to describe this solution we introduce 
some special notations. It is easy to see that at C\ = we can rewrite kinematic 
bivector ( 13. ip in the following form 

x 2 d^k A 



A = I -73 7i I (3-4) 
where antisymmetric matrix <9 7 A is defined by 



—A <9 7 A 
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t^r-t^,-^) (3.5) 
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In the similar notations at C\ — the second bivector is equal to 

A' 





( x 2 <9 7 A 
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where 



2x 2 (/3 1 M 1 2 + /3 2 Mf) , , 

a = -— -= , ft = x di-1, = x a 2 - 1 . 

(A - /3 2 )7 3 

Matrix IT is equal to 

2 X 2 / /3iAfi(M 2 73-72M 3 ) -^i(72M 2 M 3 +73M 1 2 ) Pi^iM? +M|) 

n = 7^ o \ 2 ft!(7lMiM 3 +7 3 Mf) -j82M 2 (Afi73-7iAf3) -/3 2 7i(M 1 2 +M|) 

IPl - P2J73 \ /3 2 7 3 M 2 M 3 -/3 l73 MiM 3 (Pi-fchzMi.Mi , 
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and antisymmetric matrix d^H reads as 



/ M 3 M 2 \ 2><2 / /3 273 - ft72 

= K Z -M 3 Mi + — — -/3 2 73 /3 2 7i 

V -M 2 -Mi / (A — P2J73 Y ^ l72 _p 2 ^ 1 o 

Instead of antisymmetric matrix A ( 13. 4 p in the second bivector we have symmetric 
matrix 

f -73 72 \ 
A ' _ -73 71 

27172 
72 7i , 

V 73 / 

whereas definition of <9 7 A' is completely similar to (13. 5p 

At x — > one get bi-hamiltonian structure for the Chaplygin system, which has 
been obtained in [IB]. We believe that bivector P' has some algebro-geometric 
justification, similar to compatible bivectors on so(n) from [2]. 

Remark 2 Usually the second Poisson bivector P' is the Lie derivative of P 
along some polynomial Liouville vector field X 

P' = £x{P), 

see [HI [131 [HI [16] . F° r the Bogoyavlensky system we could not find such Liouville 
vector field. So, we can not say that bivector P' (13.61) is the 2-coboundary asso- 
ciated with the Liouville vector field X in the Poisson-Lichnerowicz cohomology 
defined by P. 

To sum up, using applicable ansatze for the Liouville vector field X we get a 
real relatively simple quadratic bivector ( 13. 6 p and some more complicated complex 
bivectors. Modern computer software allows to do it on a personal computer 
wasting only few seconds. The application of this Poisson bivector will be given 
in the next section. 



4 Variables of separation and separation rela- 
tions 



The second step in the bi-hamiltonian method of separation of variables is calcu- 
lation of canonical variables of separation . . . , q n ,pi, . . . ,p n ) and separation 
relations of the form 



i(qi,Pi,Hi, . . .,H n ) = , i = 1, 



with det 



dHi 



^0. (4.1) 
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The reason for this definition is that the stationary Hamilton- Jacobi equations 
for the Hamiltonians Hi can be collectively solved by the additively separated 
complete integral 



W(qi, . . . ,q n ;ax, . . . ,a n ) =) j Wi{qi\ ax,..., a n ) , (4.2) 



i=i 




where Wi are found by quadratures as solutions of ordinary differential equations. 

According to [5j |13j HE] , separated coordinates qj are the eigenvalues of the 
control matrix F defined by 

P'dH = P(FdH). 

Its eigenvalues coincide with the Darboux-Nijenhuis coordinates (eigenvalues of 
the recursion operator) on the corresponding symplectic leaves. Using control 
matrix F we can avoid the procedure of restriction of the bivectors P and P' on 
symplectic leaves, that is a necessary intermediate calculation for the construction 
of the recursion operator [3]. 

In our case for the Poisson bivector P' (13.61) control matrix F reads as 



(4.3) 



where 

u = -(PiMl + f3 2 M 2 2 ) , v = -*T 2 (/3i - f3 2 )pl . 
The eigenvalues of this matrix F are the required variables of separation q 12 

qi = , Q2 = • (4.4) 

V V 

These variables has been introduced in [1] without any explanations and reason- 
able arguments. We reproduce this result in framework of the generic method 
based on direct solution of the equations (13. 2p . According to [5J, eigenvectors of 
the control matrix F form the Stackel matrix S 

F = S ( q ' Mr 1 
\ q 2 J 

whose entries S^- depend only on a pair (qi,Pi) of the canonical variables of 
separation. In our case matrix S is equal to 



S 



1 1 

1\J Pi 2 — Pt 2 



It means that we have non-Stackel integrable system with non-affine in Hi or H 2 
separated relations (14. ip . similar to the generalized Chaplygin system [16J and 
the Kowalevski top [T7] . 



6 



From the definitions of separation coordinates gi )2 (I4.4p and H 2 (12. 5p we 
immediately obtain 



V 2/3 1 (g 1 -g 2 ) 2 V 2/5 2 (g 2 - gi ) 2 

and 



73 



2x# 2 1/4 



v^A-AOfo-gi)' 

Such as Ci = and 

4/^^(1 + g fe ) 4/3 27l M 2 (l - g fe ) 



g fc = {H, q k } 



we have 



7i 



72 



73 73 

(1 + g 2 ) gi - (1 + gi) g 2 

4(gi-g 2 )v / /32(A-/3 2 ) v / ( 1 + ?i)( 1 + ?2) ' 

X (1 -92)9l - (1 ~gl)g2 

4(gi-g 2 )v / /3i(/3 2 -/3i) v /(f-g 1 )(f-g 2 ) ' 

g!(f-g 2 ) + g 2 (l-g 2 ) 
4(gi - g 2 )>/-/3i/32(l-?i 2 )(l-92 2 ) ' 



Substituting these expressions into the Hamiltonian H\ and Cazimir C 2 and solv- 
ing the resulting equations with respect to gi )2 one gets a pair of the Abel-Jacobi 
equations 

X ^ ' 1. (4.5) 



^(^-^(Aigi+z/i) ' v/4(g 2 2 _ i)(A 2 g 2 + /i2 ) 
where 

Ai j2 = (/3i-/3 2 )(^2(/3i + /32) + x 2 (/J 1 ± v ^ 2 )), 

/ii, 2 = (/3i-/3 2 ) 2 C 2 + x 2 (/3 1 +/3 2 )(F 1 ± v 7 ^)- 

Let us note that ( 14. 5ft are degenerate Abel-Jacobi equations, i. e. each of them 
depends on a unique variable qi or g 2 only, and a two-dimensional Abel torus 
splits into one-dimensional tori. 

According to [9], the remaining separation variables pi ;2 are equal to 

Pk = —. 2W , q \ . ry, {quPk} = S ik , i,k = 1,2. 

4 (! - %) (/M 1 + Ik) + 02{1 ~ qk)) 

They satisfy to the following separated relations which directly follow from the 
Abel-Jacobi equations (j4.5|) . 

$i, 2 (g,p) = 4x 2 (g 2 - 1) (ft(l + g) + 2 {1 - g)) V - A x , 2 g - /ii, 2 = , (4.6) 



here q = gi )2 and p = pi j2 . At x — > these equations coincide with the separated 
relations for the Chaplygin system, see [T6] . 

The third part of the Jacobi method consists of the construction of new 
integrable systems starting with known variables of separation and some other 
separated relations. Namely, if we substitute our variables of separation g 12 and 
Pi y2 into the following deformation of (14.61) 

<5> {d) {q,p) = ® 1 $ 2 -d l q-d 2 = 0, d 1} d 2 eR, (4.7) 

and solve the resulting equations with respect to integrals of motion #1,2, then 
we get rational generalization of the initial polynomial Hamilton function 



.2 



- ftx 2 (<i, + d 2 )(ft - ft)A/ 2 2 + (/3f(d! - d 2 ) + + <i 2 ))7 s 2 
If di = d 2 ((3i — fi 2 )/{fii + (3 2 ) this Hamiltonian looks like 

rr( d ) rr 2d 2 H 2 (3^(3, - (3 2 ) 

n 1 — tlx + 



(/3i + /3 2 )(x 2 M 1 2 + x 2 M|+ 7 3 2 ) 



At x — ^ we obtain the Hamilton function for the generalized Chaplygin system 
studied in [TB]. The main problem of this part of the Jacobi method is how to 
get the Hamiltonian to be interesting to physics. 



5 Conclusion 

Starting with the integrals of motion for the Bogoyavlensky system on so(4) we 
found polynomial in momenta Poisson bivector P', which are compatible with the 
canonical Poisson bivector P on zero-level of the Casimir function C\. Then in 
framework of the bi-hamiltonian geometry we reproduce known separation vari- 
ables and separated relations. Some rational generalization of the Bogoyavlensky 
system is considered. 

This example may be useful for creating a general theory, which takes the 
constructive answers on the main open questions: 

• how to get the Poisson bivectors P' on so(n) compatible with P ; 

• how to describe all the natural Hamilton functions associated with a given 
P'. 

Now we have some particular answers obtained by direct tedious computations 
only [21[ni [13]. 

The author wish to thank A.V. Tsiganov for formulation of the problem and 
stimulating discussions. 
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